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LINE-OF-SIGHT VELOCITY AS AN INDICATOR OF PERICYNTHION 
ALTITUDE OF SELENOCENTRIC ORBITS 

By James Lloyd Williams 
Langley Research Center 

SUMMARY 

An investigation was made to determine the feasibility of using simplified equations 
to predict the pericynthion of a selenocentric orbit. It was assumed that the angular rate 
of the line of sight with respect to a prominent lunar surface feature was available for 
positions along a ballistic trajectory. The equations used were derived in terms of this 
variable by utilizing the linearized rendezvous equations and certain parameters of a 
circular selenocentric orbit. The procedure used in the evaluation of these equations 
was to calculate the pericynthion of selected reference Hohmann trajectories and to com- 
pare the values calculated by these equations with the known reference orbit values. 

The results of the computations indicate that 

(1) The pericynthion altitude predicted by the approximate equations generally is 
not accurate. 

(2) The predicted pericynthion altitude is very sensitive to the reference orbit 
altitude, to the times at which measurements of the rotation of the line of sight are made, 
and to the length of time between the measurements of the angular rate of line of sight. 

INTRODUCTION 

A critical phase of the lunar mission is the establishment of a safe elliptical trans- 
fer orbit between the selenocentric circular orbit of the command module and a pre- 
scribed position close to the lunar surface, the pericynthion of the transfer orbit. The 
final touchdown maneuver is initiated from the pericynthion of the transfer orbit. Moder- 
ate errors in the transfer maneuver could result in an unsafe orbit (impact or near impact 
trajectory) in this phase of the lunar mission. It is apparent that mission reliability could 
be enhanced if simplified methods of determining orbit pericynthion were available in the 
event of automatic guidance failure. 

Several approaches to simplifying orbit determination, with varying degrees of 
complexity, have been reported in references 1 and 2. In reference 1 equations based on 
various combinations of orbital parameters are presented, and in reference 2 errors in 



pericynthion based on altitude -deviation measurements are evaluated. The approach 
studied in this investigation made use of the line-of-sight rotational velocity in respect . 
to a prominent lunar surface feature for the determination of orbit pericynthion. 

This report presents linearized equations for pericynthion prediction which were 
derived in terms of the line-of-sight rotational velocity, relative to lunar surface fea- 
tures, by making use of linearized rendezvous equations. Pericynthion altitudes calcu- 
lated from these equations for three known selenocentric elliptic orbits are also pre- 
sented. The trajectories selected for these calculations had pericynthions of -186, 

10 862, and 27 611 meters. 


SYMBOLS 

h altitude above lunar surface, meters 

Ahp numerical difference between true pericynthion and predicted pericynthion, 

meters 

m mass of body 

r radial distance from center of lunar sphere to vehicle, meters 

r radial velocity component of vehicle, meters/second 

r m radius of lunar sphere, meters 

rd velocity component of vehicle perpendicular to radius vector, meters/second 

t time measured from apocynthion, seconds 

x s relative displacement measured along circular selenocentric orbit from 

origin of moving axis, meters 

y s relative displacement measured along radius r from end of x s displace- 

ment to landing vehicle, meters 

x relative rectangular displacement perpendicular to radius r 0 and measured 

from origin of moving-axis system, meters 

y relative rectangular displacement measured along radius r Q from origin of 

moving-axis system, meters 
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X,Y,Z moving rectangular coordinate system 

6 ' angular position of vehicle measured from reference position, degree or 

radian 

X angular separation of moving-axis -system origin and vehicle in orbit 

p distance from landing vehicle to prominent lunar surface feature, meters 

cu angular rate of moving-axis origin in selenocentric orbit, radian/second 

co^g angular rate of line of sight, radians/second 

<p angle between vehicle's local vertical and line of sight to surface feature, 

degrees 


Subscripts: 

p pericynthion 

1,2,3 order of reading 

o selenocentric circular orbit conditions 

Dots over symbols denote derivatives with respect to time. 

ANALYSIS 

It is desired to ascertain whether orbital characteristics can be obtained from the 
angular rate of the line of sight with respect to a prominent lunar surface feature. The 
primary coordinate system used in the derivation of the equations was a cylindrical shell 
coordinate system the origin of which was assumed to be in a circular selenocentric orbit. 
A sketch of this system denoted by the displacement components x s ,y s can be seen in 
figure 1, along with a typical moving rectangular coordinate system X,Y. A set of 
transformation equations that connects the shell and rectangular coordinate systems is 
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Figure 1.- Sketch of coordinate axis systems, angles, and displacements. 



A typical representation of the geometrical relationships between the line -of -sight and 
velocity components is shown in figure 2. The vehicle velocity normal to the line of 
sight can be expressed as follows: 

u>l s p = (r0)cos (p + r sin 0 (3) 

The velocity components, in polar form, can be expressed as 

(r 6) = x cos A - y sin A + r 0 u> cos A 


r = x sin A + y cos A + r Q oj sin A 
Equation (3) therefore can be written as 

o>£ S p = (x cos A - y sin A + r 0 w cos a)cos (p + (x sin A + y cos A + r 0 u> sin A^sin 4> (4) 

Differentiation of equations (1) and (2) yields 

x = y s sin ^ + (r Q + y s )(cos 
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Figure 2.- Sketch showing angular displacements and angular and linear rates. 


y = y s cosp£-(r ot y s )^in|a^£ 


Use of these relations in equation (4) results in 


w zs p = ys sin + ( r o + Ys) ^ cos cos x cos <P 


- y s cos ^ - (r 0 + y s )^ sin sin x cos 

+ r 0 co(cos X cos cp + sin X sin 0) + y s cos ^ - (r 0 + y s )|sin C os X sin <p 


. . Xc 

+ y s sin -a + 
r o 


( r o + y^ cos sin * sin 0 


Use of the relation X = -§• in equation (5) results in 

r o 


co£ s P = y s sin (p + ^r Q + y s ^^ cos <p + r n u>(cos X cos (p + sin X sin 


It is assumed that cp and X are small angles so that sin <p = <p, cos <p = 1; sin X = X, 
cos X = 1. For this investigation the maximum separation between the origin of the 
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moving axis and the landing vehicle X was about 6°. This assumption implies that 
p ~ h Q + y s . Under these assumptions, equation (6) becomes 

w zs( h o + y s ) = y s 0 + ( r o + y s )?J + r o w (* + x 4>) ( 7 ) 

If it is assumed that r Q + y s ~ r Q and second-order terms are neglected, equation (7) 
becomes 




+ r o w 


( 8 ) 


Based on linearized equations of motion similar to the method developed in refer- 
ence 3, the position and velocity of a mass in a nonthrusting trajectory relative to a 
moving-axis system in a circular orbit are given by the following equations: 


x s(t) = 


4x s (0) 


CO 


+ 6yo(0))sin cot + 


2y s (0) 


CO 


cos cot - (6y s (0) + 


3x s (0)' 


CO 


cot + x s (0) - 


2y s (0) 


CO 


(9a) 


= + 6 y s (°)) cos wt - sin - 6 y s (°) - 


(9b) 


y s W = 


2x s (0) 


CO 


+ 3y g (0)Jcos cot + sin cot + 4y s (0) + ■ 


(9c) 


YsW 

CO 


2x s (0) 


CO 


3y s (0) sin wt 


y s (°) 


CO 


cos cot 


(9d) 


where x s (0), y s (0), x s (0), and y s (0) refer to conditions at specified positions along 
the ballistic trajectory. These equations differ in sign from those in reference 3 because 
of the clockwise rotation of the moving-axis system used in this investigation. With the 
aid of equations (9b) and (9c), equation (8) can be written as 


A(t) 


x s (°) 

co 


+ B(t)y s (0) + C(t) 


y s (°) 


= D(t) 


(10) 


where A (t), B(t), C(t), and D(t) are defined as follows: 



An equation that expresses the maximum displacement in the y s -direction can be 
obtained from equations (9c) and (9d). This condition implies y s (t) = 0 and from equa- 
tion (9d) the following equation is obtained 


-y s (o) 


cet = tan 


-1 


OJ 


i2x s (0) 


(ID 




+ 3y s (0)/ 


The substitution of equation (11) for cot into equation (9c) results in the following equa- 
tion for the maximum displacement in the negative y s (t) direction denoted by (y s ) p 


( y s) P = + 3 y«(°)) cos 


tan 


-1 


I :m 


2x s (0) „ . x 

s + 3y s (0)i 


't 


co 


y s (0) . 

+ — sin 


co 


/ -y s (»> 


tan“lj 


co 


2x s (0) 


+ 3y s (0)i 


2x s (0) 


+ 4y s (0) + — | 


( 12 ) 


The pericynthion of the orbit is 


h P “ ( y s)p + h o 

Equations (10) to (13) are used to predict orbit pericynthion based on values of angular 
velocity of the line of sight with respect to a lunar surface feature. 

To determine the three variables x s (0), y s (0), and y s (0) in equation (10) 
requires that three consecutive values of < j ^ i s and associated times t in orbit be 
available. The magnitude of coj s used in the derived equations was calculated from the 
expression coj s = r0/h (for 4 > = 0) at selected times along the known elliptical transfer 
orbits. The times were measured from the apocynthion of the orbits. The coeffi- 
cients A, B, C, and D were evaluated, after which standard procedures for solving 



simultaneous equations were used to obtain the relative velocity and displacement com- 
ponents x s (0), y s (0), and y s (0). These values were assumed to be the conditions at the 
time of the first o^ s measurement and were used to predict pericynthion altitude hp 
from equation (13) with the aid of equations (11) and (12). 

RESULTS AND DISCUSSION 

Three ballistic trajectories were selected for use in the evaluation of the equations 
derived herein. Each had an apocynthion of 148 160 meters and had pericynthion altitudes 
of 27 611 meters, 10 862 meters, and -186 meters. Values of u>£ S were calculated along 
each of the trajectories (fig. 3) and these were used in subsequent computations. Since a 
number of approximations and assumptions were made in the derivation of the equations, 
it was realized that there would be errors involved in the predicted pericynthion altitudes. 
During this study the effects of the following variables were examined: time of assumed 
ci>£ s measurements, reference orbit altitude, and error in estimated initial conditions. 



Figure 3.- Variation of «j s for three selected descent ballistic trajectories. 0-0. 

Effect of Assumed Time of u )^ s Measurements 

The reference circular selenocentric orbit, which contains the origin of the moving- 
axis system, passes through the apocynthion of the transfer orbit. Time was assumed to 
be measured from apocynthion. The computed initial conditions corresponded to the con- 
ditions at various times from apocynthion. As mentioned previously, these initial values 
were assumed to be the conditions at the first a > is measurement. The predicted peri- 
cynthion altitudes for the three reference descent orbits are shown in figure 4 plotted 
against time measured from apocynthion to the first co^ g value. The specific times 
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used for the assumed <x>i s measurements are given in table I. These times were 
chosen arbitrarily to extend from apocynthion to near pericynthion (about one-half the 
orbital period). A study of figure 4 shows that the trends are the same for all three 
orbits, and generally indicate pericynthion altitudes lower than the reference -orbit peri 
cynthions. However, these results indicate an improved prediction as the time for mea 
surement is delayed. 



Figure 4.- Predicted pericynthion altitude plotted against time of first value. 
Reference circular orbit altitude of 148 160 meters. 
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TABLE I.- PREDICTED hp CALCULATED FROM THREE VALUES OF w ls 
AND THE ERROR Ah p AS A FUNCTION OF VARIOUS ORBITAL TIMES 
[Reference circular orbit altitude of 148160 meter^j 



t,seconds,measured from apocynthion 



hp, reference, m 




hp, predicted, m 

Ahp, m 

*i 

*2 

'3 





383 

701 

1 274 

26136 

-1475 


701 

1 108 

1718 

22752 

-4859 


1 108 

1718 

2421 

22440 

-5171 

2761 1 

1575 

1859 

2296 

22835 

-4776 


1859 

2296 

2655 

25908 

-1703 


2145 

2533 

3193 

28477 

866 


2315 

2655 

3193 

28680 

1 069 


361 

863 

1459 

6079 

-4783 


658 

1 186 

171 1 

3028 

-7834 


863 

1 186 

1 586 

2245 

-8617 


1 186 

171 1 

2551 

3348 

-7514 

10862 

1459 

171 1 

2182 

2056 

-8806 


171 1 

2098 

2452 

5112 

- 5750 


2082 

2452 

2921 

9444 

- 14 18 


2358 

2659 

3106 

10893 

31 


345 

826 

1 391 

-6488 

-6302 


629 

1 1 34 

1854 

-II 301 

-1 1 1 1 5 


989 

1391 

1854 

- 1 2 520 

-12334 


1 1 34 

1 854 

2465 

-9654 

-9468 

-186 

1511 

1968 

2377 

-9348 

-9162 


1627 

2058 

2465 

-8099 

-7913 


1854 

2213 

2658 

-5508 

-5322 


1982 

2465 

2898 

-3049 

-2863 


2465 

2658 

2898 

-1 4 28 

-1242 


Some calculations were also made for a pericynthion altitude of 10 862 meters in 
which the time interval between a>i s measurements were approximately constant. The 
results of these calculations are presented in figure 5, as predicted hp values plotted 
against time of the first measurement. Three approximately constant time intervals 

of 200, 600, and 1000 seconds are shown in figure 5. For purposes of comparison, the 
results for the arbitrary intervals (nonconstant) are also presented. The nonconstant time 
intervals were generally in the vicinity of 200 to 600 seconds. A study of this figure shows 
that predicted hp values closer to the true value of hp were obtained by lengthening the 
time interval between assumed a^ g measurements. (Compare results for 200, 600, and 
1000 seconds.) It should be noted that this effect occurred primarily at the larger values 
of delay time from apocynthion. This trend suggests that predicted hp values in the 
neighborhood of the true values can be obtained by a proper selection of measure- 

ment as a function of delay time from apocynthion. However, it should be pointed out 
that the total time necessary to obtain the necessary three c o^ s measurements for a 
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reasonably accurate prediction 
required that the last measured 
value occur so close to the time of 
pericynthion that sufficient time for 
corrective measures may not be 
available in the event that a danger- 
ously low altitude was predicted. 

The error in predicting peri- 
cynthion altitude, based on the 
results in table I, is shown in fig- 
ure 6. This figure presents Ahp 
plotted against time measured from 
apocynthion to the first c o^ s value. 
The error in predicted pericynthion 
altitudes varies inversely with actual 
pericynthion altitude. This variation appeared reasonable since the linearized rendezvous 
equations have an inherent assumption that gravity varies linearly from the reference 
orbit altitude. This result seems to indicate that a reference altitude closer to the aver- 
age orbit altitude would be better than one at either extreme of the actual orbit. The mag- 
nitude of the error can be reduced somewhat by lengthening the intervals between co^ g 
measurements. This result can be inferred, for 10 862-meter pericynthion orbit, from a 
comparison of the results in figures 5 and 6. 


Figure 5.- Variation of predicted pericynthion attitude with time plotted 
at first Wj s value for several constant time intervals between 
assumed measurements. Reference selenocentric orbit, 148 160 meters. 



Figure 6.- Error in predicted pericynthion altitude plotted against time of first 
value. Reference circular orbit altitude of 148 160 meters. 
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Effect of Reference Orbit Altitude 

A series of computations were made with the origin of the moving-axis system in a 
circular orbit at an altitude of 74 080 meters, that is, one-half of the maximum orbit alti- 
tude. It was further assumed that the landing vehicle initiated its descent maneuver from 
the maximum orbit altitude of 148160 meters. The results are shown in figure 7, which 
shows the predicted pericynthion altitudes plotted against time measured from apocyn- 
thion to the first o>£ S value. 



Figure 7.- Predicted pericynthion altitude plotted against time of first wj s value. 
Reference circular orbit altitude of 74 080 meters. 
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Table n indicates the times at which the co^ s values were obtained. A comparison of 
figure 7 with figure 4 shows that decreasing the altitude of the reference selenocentric 
orbit did not appreciably affect the trends of the variation of h p with time of measure- 
ment; it did, however, shift the entire pattern upward toward higher values of h p . Thus, 
the predicted values of h p were higher than the true values at small delay times and 
somewhat closer to the actual values at the larger delay times than the results for the 
higher reference orbit. 

Results for calculations made with approximately constant time intervals between 
readings are presented in figure 8 for a pericynthion altitude of 10 860 meters. This fig- 
ure presents predicted h p values plotted against time measured from apocynthion for 
constant time intervals of about 200, 600, and 1000 seconds. Also included in this figure 
are the results for nonconstant time intervals. A comparison of figure 8 with figure 5 


TABLE D.- PREDICTED h p CALCULATED FROM THREE VALUES OF w Zs 
AND THE ERROR Ah p AS A FUNCTION OF VARIOUS ORBITAL TIMES 
[Reference circular orbit altitude of 74 080 meters] 


hp, reference, m 

t, seconds,measured from apocynlhion 

hp, predicted, m 

Ah Pl m 


*1 

*2 

*3 



383 

701 

1274 

36170 

8559 


701 

1108 

1718 

30280 

2669 


1 108 

1718 

2421 

25536 

-2075 

2761 1 

1575 

1859 

2296 

25145 

-2466 


1859 

2296 

2655 

26224 

-1 387 


2145 

2533 

3193 

2 7 745 

134 


2315 

2655 

3193 

27868 

257 


361 

863 

1459 

17671 

6809 


658 

1186 

1711 

12483 

1 621 


863 

1186 

1586 

1 1 748 

886 

10 862 

1 186 

1711 

2551 

7078 

-3784 

1459 

171 1 

2182 

6543 

-4319 


171 1 

2098 

2452 

7264 

-3598 


2082 

2452 

2921 

9690 

-1172 


2358 

2659 

3106 

10710 

-152 


345 

826 

1 391 

6968 

7 154 


629 

1 134 

1854 

-754 

-568 


989 

1391 

1854 

-3631 

-3445 

-186 

1 1 34 

1854 

2465 

-5159 

-4973 


1511 

1968 

2377 

-5625 

-5439 


1627 

2058 

2465 

-5147 

-4961 


1854 

2213 

2658 

-3813 

-3627 


1982 

2465 

2898 

-2291 

-2105 


2465 

2658 

2898 

-1132 

-946 
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again shows the upward shift and the similarity of the trends of variation with time of 
measurement. These characteristics have been discussed previously for figure 5. The 
data of figure 9 show the error in predicted h p altitude as a function of the time of the 
first measurement. The effects of decreasing the orbit altitude to 74 080 meters 

(which were pointed out previously) are shown in this figure. As pointed out earlier, 
some reduction in the magnitude of the error can be obtained by lengthening the time 
interval between co^ g measurements. (Compare results in figs. 8 and 9 for the 
10 862-meter pericynthion.) 



Figure 8.- Variation of predicted pericynthion altitude with time 
plotted at first value for several constant time intervals 
between assumed measurements. Reference selenocentric 
orbit, 74 084 meters. 



Figure 9.- Error in predicted pericynthion altitude plotted against time of first 
Wj s value. Reference circular orbit altitude of 74 080 meters. 
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Sources of Errors in Predicted Pericynthions 

Errors in predicted pericynthion were very sizable in some cases. There appeared 
to be two possible major sources of error: the estimated initial conditions, and the use of 
the approximate equations for estimating pericynthion altitude. 

Figure 10 compares the exact and the estimated initial conditions (based on calcu- 
lated u> ig 

and estimated values for the initial conditions diverge as tj increases, and suggest the 
increased importance of second-order terms neglected in the derivation of equation (10). 
The divergence in estimated and exact initial values varies regularly and, therefore, 
would not account for the type of Ahp variation shown in figures 6 and 9. 

At this point the exact initial conditions (values obtained from exact orbit equations) 
were used with the pericynthion prediction equations (that is, equations (11), (12), and (13)), 
and results are shown in figure 11. This figure shows a comparison of estimated hp 
obtained by using these exact values and by using estimated initial conditions based on 
o^ g measurements. It is apparent that use of the exact initial conditions in the approxi- 
mate equations generally does not improve the estimated pericynthion prediction. To 
some extent the use of exact initial conditions made matters worse in that there is no 
longer a convergence toward the proper value as the time for co^ g measurements is 


values) at the time of the first reading. The differences in these exact 



Figure 11.- Comparison of pericynthion altitudes calculated from exact and estimated 
initial conditions for orbit with 10 862-meter pericynthion. Linearized rendezvous 
equations were used for exact calculations. Reference circular orbit altitude was 
148 160 meters. 
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delayed; the value to which the predictions appear to converge overestimates pericynthion 
altitude. These trends suggest that perhaps better results could be obtained with this pro- 
cedure measurement) if an improved set of linearized equations could be developed. 

CONCLUDING REMARKS 

An investigation was made to determine the feasibility of using simplified equations 
to predict orbit pericynthion. It was assumed that the angular rate of the line of sight 
with respect to a prominent lunar surface feature was available for positions along seleno- 
centric orbits. The equations used were derived in terms of this variable by utilizing the 
linearized equations of relative motion of orbiting bodies and certain parameters of a cir- 
cular selenocentric orbit. The procedure used for the evaluation of these equations was 
to calculate pericynthion altitudes of known reference orbits and to compare the estimated 
values with the actual values. 

The results of the computations indicate that 

(1) The pericynthion altitude predicted by the approximate equations generally is 
not accurate. 

(2) The predicted pericynthion altitude is very sensitive to the reference orbit alti- 
tude and to the times (location in the transfer orbit) at which measurements of the rotation 
of the line of sight are made, and to the length of time between the measurements of the 
angular rate of line of sight. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 31, 1967, 

125-17-05-01-23. 


REFERENCES 

1. Gersten, Robert H.; and Schwarzbein, Z. E.: Self-Contained Orbit Determination 

Techniques. [Paper] No. 63-431, Am. Inst. Aeron. Astronaut., Aug. 1963. 

2. Reid, Richard: A Simplified Technique for Determining Deviation in the Lunar 

Transfer Orbit Ephemeris. NASA TN D-1837, 1964. 

3. Eggleston, John M.; and Beck, Harold D.: A Study of the Positions and Velocities of a 

Space Station and a Ferry Vehicle During Rendezvous and Return. NASA TR R-87, 
1961. 


NASA- Langley, 1968 30 L-5296 


17 



